NON-FRIEDRICHS SELF-ADJOINT EXTENSIONS OF THE 

LAPLACIAN IN R d 

PAUL LIN 

Abstract. This is an expository paper about self- adjoint extensions of the 
Laplacian on R d , initially defined on functions supported away from a point. 
Let L be the Laplacian — Xh=i 7TT w ith domain C^°(IR d \{0}). We determine 



(N 



1 dx\ 

all the self-adjoint extensions of L by calculating the deficiency subspaces of 
the closure L of L. We give the details of the calculations here since, to my 
knowledge, they are nowhere else in the literature. In the case d = 3, the self- 
adjoint extensions Lg are parametrized by a circle 9 € [— tt,tv), or equivalently 

■ by the real line plus a point fi E KU {oo}. We show that the non-Friedrichs 
extensions have either a single eigenvalue or a single resonance. This is of some 
interest since self-adjoint Schrodinger operators A + V on R d either have no 

, resonance (if V = 0) or infinitely many (if V ^ 0), see RBM ( Prop. 4.3]. 

■ For these operators , the kernel of (L M — A 2 ) -1 is calculated explicitly, see 
[AGHH . Finally we use this knowledge to study a wave equation involving L^, 

^3 ■ and explicitly determine the wave kernel. 



1. The operator L and its closure 



> 

We work in the Hilbert space L 2 (R d ), and start with the unbounded operator 
L = -J2t=i in with domain D(L) = C~(M d \{0}). Note that L is symmetric, 



1 dx i 

which can be shown by integration by parts, but it is not self-adjoint as the domain 



q ! of its adjoint D(L*) contains C^°(R d ) which is strictly bigger than D(L). Recall 
that the deficiency subspaces Jf + , JPL of L are the null spaces of the operators 
i — L*, i + L* respectively. By applying von Neumann's theorem to L, we know 
that the two deficiency subspaces of L have the same dimension hence L has self- 
^ \ adjoint extensions. To find out how many self-adjoint extensions L has and what 
they look like, we use the following proposition due to von Neumann, for the proof 
see (ESI Sec. X.l]: 

Proposition 1.1. Let A be a closed symmetric operator with equal deficiency 
indices. Then there is a one- one correspondence between self- adjoint extensions 
of A and unitary maps from Jtf + onto J^- . If U is such a unitary map, the 
corresponding self-adjoint extension Ajj has domain 

D(Au) = & + ip + Utjj : <p G D(A),ip e 

and 

Au((p + tp + Uip) = Acp + itp- iUip. 



Key words and phrases. Laplacian, non-Friedrichs self-adjoint extension, eigenvalue, 
resonance. 
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The above proposition requires a closed symmetric operator. To find the domain 
of L means to complete C^°(lR a! \{0}) under the norm 1 1 • 1 1 + 1 |L(-)| |. (In this paper, 
the norm notation || • || always denotes the L 2 -norm.) This norm is equivalent to 
the W 2 > 2 -norm hence we have D{L) = W 2,2 (R d \{0}) C W 2 > 2 (M. d ). The situation 
differs as dimension varies. We first treat the case d > 4. 

1.1. Case d>4. 

Proposition 1.2. ForJBL d with d > 4 ; we have 

D{L) = W 2 ' 2 (R d ). 

Proof. From above, we already know that D(L) C W 2,2 (M. d ). Since C^°(lR d ) is 
dense in W 2,2 (M. d ), we pick any if> e C^°(M. d ), and approximate it with a sequence 
of functions in C^°(lR d \{0}) that converges to ip under the W 2 ' 2 -norm. Let (p be a 
smooth function such that (p — 1 on Bi(0) and ip — outside -Bi(O). Then define 
ip € (x) = 1 — <£>(-). We now show that ^<£> e converges to if) under the W 2 ' 2 -norm 
when e approaches 0. Indeed, ip(p € £ (M, d \{0}) , and 

||^ e -V;|| 2 = / \iP(x) V( -)\ 2 dx < Ml [ \ V {-)\ 2 dx = e d M 2 |M| 2 ->• 0, as e -> 0, 

where Mi is the maximum value of ■0. Now consider \\L(if)ip e — if)))\\, product rule 
creates three terms, 

||(^)(v? e -l)|| 2 = f \(Li>)(xM-)\ 2 dx < M% [ \ V (-)\ 2 dx = e d Ml\\ V \\ 2 ^0, as e 

JR d 6 JR d e 

where M 2 is the maximum value of Lip. Since d — 2 > 0, 
||V^-V(^-1)|| 2 <M 2 / \V V (-)\ 2 dx = e~ 2 Ml [ \iy V )(-)\ 2 dx 

JR d e JR d e 

= e d - 2 M 2 |||Vv?||| 2 ^0ase^0, 
the partial derivatives of if) take values smaller than M 3 . At last, since d — 4 > 0, 

||^(^-1)|| 2 <M 2 / \L V (-)\ 2 dx = M 2 ^ [ \{Ltp){-)\ 2 dx 

JR d e JR d e 

= M 2 e d - 4 \\Lcp\\ 2 ^Oase^O. 

We have established in the case d > 4, 

D{L) = W 2 ' 2 (R d ). 

□ 

1.2. Case d=4. We can see for the case d < 4, the above proof fails because 
M 2 e d_4 | |L</?| | 2 doesn't converge to 0. In fact, when d < 4, the domain of L is 
smaller than W^ 2 ' 2 (M d ). For the case d = 4, we still have £>(Z) = W 2 ' 2 {R 4 ), but we 
can't use the same ip t as defined in the above proof because | \L(ip e — 1)|| = | |L<^(^)| | 
is independent of e. So to get convergence we need to define <p e in a way that breaks 
the scaling invariance of \\L((p e — 1)||. 

Proposition 1.3. ForM 4 , we have 

D(L) = W^ 2 ' 2 (R 4 ). 
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Proof. As before, we pick any arbitrary ip e C^°(M. 4 ). Let ip : [0, oo) — > [0, 1] be 
a smooth function such that such that </?([0, |]) = 1 and (p([l,oo)) = 0. Then 

define (p e (x) := 1 — We now verify that as e — > 0, -0</? e approximates -0 

under the VU 2 ' 2 -norm. Indeed, ip<p t G C^°(IR 4 \{0}), and 



\\^-n 2 = [ hK*M(-r)i 2 ^- 



We show that the above norm tends to when n becomes large. Consider any x, 
when e < \x\, we have (^-) e > 1, hence \ip(x)(p((^-) e ) | 2 — >■ pointwise. Moreover, 
\i/j(x)ip((— ) e ) | 2 is bounded above by the L 1 function |-0| 2 , so we can conclude 
\\ip(f e — -0| | — > as e — )• 0. 

Again as before, the term \\L(ifjtp t — \ produces three terms: \\(Lxjj)(ip e — 1)||, 
||V-0 • V(y? e — 1)|| and \\ipL(ip e — 1)||. Here we just show how to deal with the 
most difficult term | \ipL((p e — 1) 1 1 where both derivatives fall on the cutoff function, 
as the computations for the others are similar, but easier. We will use the full 
strength of d — 4. As before, we just need to estimate ||£(y? e — 1)||. We will first 
calculate this norm for general 1 < d < 4, then substitute the dimension d = 4, 

||L(y. £ -l)|| 2 = / \L<p((Wy)\ 2 dx. 

After computation, we obtain the expression 

= ( rf - 2 ) el -i*ry ((-r) 
w e 1 1 , 

+ e 2 - 2 ^| 2 -y((%) + e 2 -| x |-y((%). 



e e 



We will use it in the proofs for all the dimensions 1 < d < 4. In particular here 
for d — 4, it becomes 

L^((M) £ ) =2e 1 - e |x| e "V((^) e ) +e 2 - 2e |x| 2 -V((^) e ) + e 2 -|af "V ((M)^). 

Because both ip' and ip" are bounded, we just need to work with the three terms 
e 1_e |a:| e ~ 2 , e 2_2e |x| 2e ~ 2 , e 2 ~ e \x\ e ~ 2 . Because <// and p>" vanish outside the ball 
-6(0, e), we only integrate over the ball -6(0, e). The square of these terms are 
e 2 ~ 2e |x| 2e ~ 4 , e 4 ~ 4e |x| 4e ~ 4 , e 4 ~ 2e |x| 2e ~ 4 . The biggest is e 2_2e |a;| 2,E ~ 4 , so we only need 
to estimate this one. Using polar coordinates, 

/ e 2 ~^\x\ 2 ^dx = e 2 " 2e fr 2 ^dr = e 2 ~ 2e [ r^dr = - as e -> 0. 

JB(0,e) JO JO 2 

This completes the proof. □ 

Remark 1.4. Note that this convergence is only logarithmic as a function of the 
'spread' of y?((^) e ), as one would expect. 
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1.3. Case d=3. We already know the domain of L lies in W 2 ' 2 (M 3 ). By Sobolev 
Embedding Theorem we know that W 2,2 (M?) can be embedded into (IR 3 ), there- 
fore it makes sense to talk about the value of one of these functions at a single 
point, in this case the origin. Any function in D(L) is the limit of a sequence of 
continuous functions that take at the origin under the W 2,2 -noim, hence it must 
also take at the origin. Therefore we know, D{L) C {*0 £ W 2,2 (j& 3 )\tp(0) = 0}. 
In fact, they are equal. 

Proposition 1.5. In M 3 , we have 

D(L) = {4>e W 2 ' 2 (R 3 )\ip(0) = 0}. 

Proof. We pick any arbitrary ip 6 c7?°(lR 3 ) with 0(0) = 0. We use the same 
sequence as the case d — 4 to approximate ip. The calculations are similar as 
well. The difference here is we get one less power of r when we change into 
polar coordinates, but this is compensated with the Holder condition. We show 
\\ipip e - if>\\,\\(bij))((p e - 1)||, ||W • V(<p e - 1)|| -»■ as e exactly the same 
as before. The only substantially different term is \\ipL((p e — 1)||. By the Holder 
condition, we have |^(x)| 2 < C\x\ for some constant C, it follows that 

IIW^-i)|| 2 = / \i>{x)L<p((^y)\ 2 dx 

= / \^{x)L V ({^y)\ 2 dx< c [ \x\\L<p({^y)\ 2 dx. 

JB(0,e) e JB(0,e) e 

We substitute d = 3 into ((H) to obtain expression for Lip((^-) e ), 

L<p((^Y) = 6 1 -|o;r-V((^) £ )+e 2 - 2 ^i 2 -VX(^) £ )+6 2 -^r-V(( 1 f 1 ) £ ). 

It is the same as the case d — 4 except the coefficient for the first term goes down 
by 1. (Later on when we discuss the case d = 2, this term will disappear.) The 
three terms we need to work with are still e 1— e | a; | e— 2 , e 2 ~ 2e \x\ 2e ~ 2 and e 2 ~ e \x\ e ~ 2 . 
Square these terms and multiply with we get e 2_2e |x| 2e_3 , e 4_4<E |a;| 4<E_3 and 
e 4_2e |x| 2e ~ 3 . Parallel to before, we integrate them over the ball -8(0, e), and we 
check the worst term to complete the proof, 

/ e 2 - 2 '\x\ 2 ^dx = e 2 - 2 ' f r 2 "-\ 2 dr = e 2 " 2e f r^dr = - -> as e 0. 

J B(0,e) JO JO 2 

□ 

1.4. Case d=2. Here by Sobolev Embedding Theorem, we have W 2,2 (M. 2 ) C 
C7 7 (R 2 ) for any 7 < 1. Unfortunately we don't get 1, so the increase in the Holder 
exponent is not enough to compensate the loss in one power when changing into 
polar coordinates. Instead here for any ip e D(L), we have 1^(^)1 — ~ 
■0(0)1 — C|^c| (In i^y) ^ for some constant C > 0, and for small \x\, see |MTl Sec. 
4.1]. 

Proposition 1.6. ForM 2 , we have 

D{L) = {tpe W 2 ' 2 (R 2 )\ij(0) = 0}. 
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Proof. The terms | \ip(p e — ip\ \, \ \(Lif})((p e — 1)| | and \\Vtft ■ V(y e — 1)|| are estimated 
like before. Now we work on the term \\ipL((p e — 1)||. As stated above, here we 
have 

\i/)(x)\ < C\x\(\n^-)^, 
\x\ 

for some constant C > 0. Then 

\\1>L(<p e - 1)\\ 2 = / |Vi(x)2^((M y)\ 2 dx<C 2 [ \x\ 2 ln±r\Ly>((^Y)\ 2 dx. 

e JB(0,e) Fl e 

Again as before, we substitute d = 2 into (OQ) to obtain the expression for L</?((^-) e ) , 

M(-r) = e 2 - 2e ix| 2 -v(( M ) e ) +e 2 -^r-y(( M ) e )- 

e e e 

Note that in this dimension the first term in ([T]), e 1_<E |x| <E ~ 2 v9 / ((^-) e ), is gone, which 
is good news as it was previously the biggest term. Because ip' and ip" are bounded, 
we estimate e 2 ~ 2e \x\ 2e ~ 2 and e 2 ~ e \x\ e ~ 2 . We square them and then multiply with 
|x| 2 ln4| to obtain e 4 ~ 4e |x| 4e ~ 2 In A and e A ~ 2e \x\ 2e ~ 2 In A. The second term is 
bigger, 

/ e 4 - 2e |x| 2e - 2 ln i- = -e 4 " 2e [' r^hirdr = + -. 

^B(0,e) M Jo 2 4 

This expression goes to as e approaches 0. This completes the proof. □ 

1.5. Case d=l. In this case the power of r we obtained from changing into polar 
coordinates no longer exists, instead we have Holder condition for the derivatives 
as compensation. This means the domain is different from the above cases as we 
must impose condition on the derivatives. Here, by Sobolev Embedding Theorem, 
not only W 2 > 2 (R) C C§(R) we also have W l > 2 {R) C C*{R), which means it makes 
sense to talk about the derivative of a function in the domain at a single point 
as well. We will skip the proof of the following proposition as it is similar to the 
other cases. 

Proposition 1.7. For IR, we have 

D(L) = E W 2 ' 2 (R)\ip(0) = 0,^(0) = 0}. 

2. Self-adjoint extensions of L 

After finding the closures we are now in the position to calculate the self-adjoint 
extensions. 

2.1. Case d>4. In this case L is already self-adjoint. We know that by calculating 
the deficiency indices of L. Suppose that u E it means, 

(it, (L + i)f) = for all / E W 2 ' 2 {R d ). 

In Fourier space, we have (u, (|£| 2 + i)f) = for all / E W 2 > 2 (R d ). Then ((|£| 2 - 
i)u, /) = for all / E W 2 ' 2 (R d ). Since ^(W 2 ' 2 (R d )) is dense in L 2 (R d ), we have 
(|£| 2 - i)u = 0. It follows that u = 0, ie Jf+ = {0}. Similarly, Jf_ = {0}. 
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2.2. Case d=2, 3. Here for u G we have 

(u, (L + i)f) = for all / G iy 2,2 (M d ) such that /(0) = 0. 
Again we work in Fourier space, then the condition becomes 

(u, (|£| 2 + i)f) = for all / G (1 + |C| 2 ) -1 L 2 (M d ) such that / = 0. 

Let (y? be a fixed function in W 2,2 (M. d ) such that L d <p(£)d£ = 1. Then for any 
/ G H/ 2 ' 2 (M d ), the function /-/ Rd f(x)dx0 is in (l + |^| 2 )- 1 L 2 (M d ), and its integral 
over M. d is 0. Hence we have 

U(\£\ 2 + i)(f- [ f(x)dx0))=O. 

That means 

(m 2 -i)u-((\t;\ 2 -i)u,<p)J)=o. 

Since H /2 ' 2 (M d ) is dense in L 2 (IR d ), we conclude that (|£| 2 — i)u = c where c = 
(G£| 2 — ^)%0) i s a constant. Thus u = so we know is the one dimen- 

sional complex space spanned by J 7 ~ 1 ( ^l_ i )- Similarly, J£l is the one dimensional 
complex space spanned by J r ~ 1 ( i|p+| )- Since both and J£l are one dimen- 
sional, and ]|p4^ have the same norms, the unitary maps from onto 
can be parametrized by the unit circle in the complex plane, such that for each 
9 G [— 7r,7r), J r ~ 1 ( |^|2_j ) is mapped to e t6> J-'" 1 ( i|pq^ )- Therefore we have, 

Proposition 2.1. TTie self-adjoint extensions of L in W 1 , d = 2,3, can be parametrized 
by a circle 9 G [— tt, tt), with 

D(L e ) = {^ + ^-\—L_) + e w ^-\-^—) : <p G W 2 ' 2 (R d ), <p(0) = and (3 E C}, 
Iff-* Ifl +* 

and 

L,(^W- 1 (^-)+e i ^- 1 (^— )) = ). 

Remark 2.2. Notice that rast^ Tfp+i are * n L 2 (M d ) iff d < 4, hence they cannot 
be in the deficiency subspaces of L for any d > 4. 

Among all the self-adjoint extensions of L, L.^ is the special one. We denote it 
by A, and call it the Laplacian. It can be easily shown that 

D(A) = W 2 ' 2 (R d ), d = 2,3. 

It is the only self-adjoint extension whose domain is contained in the form domain 
of the closure of the quadratic form q associated with L: 

q((p,ip) = / V^p-Vipdx. 

JWL d 

Note that D(A) = W 2 ' 2 (R d ) is contained in Q(q) = H/ 1 ' 2 (R d ), it follows from [ESI 
Sec. X.3] that A is the Friedrichs extension. 
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For d = 3, let's look at another way to parametrize these self-adjoint extensions 
by considering the Taylor expansions at the origin of the functions in the domains 
of these extensions. The parameter is the ratio between the constant term and 
the coefficient of the A term in the expansions. For that, we use the following 
well known result, 

1 1 e~ A ' a '' 

(2) ^ 1( |£| 2 + A 2) = 4^^' 

for A G C\iR. From this we obtain these two expansions at the origin, 

1 1 p~cis(-f )|as| i i 

^ 1 (^2 L -) = 7~ n = TTT - T^ e_f 1 + °(M)> 

|£| 2 -r 47r |x| 47r|a;| 47r vi 17 

~ i/ 1 n 1 e~ cis ^ x \ 1 1 nf . IN 

|4 1 + « 47T \x\ 47T|X| 47T 

We will use the above expansions in the proof of the following proposition, 

Proposition 2.3. The self-adjoint extensions of L in R 3 can be parametrized by 
/iGlU {oo} with 

D(L") = W = +fi ) + ( p': ( p'e W 2 ' 2 (R 3 ), <p'(0) = 0, /3 G C}, 



where \ is a function in C^°(1R 3 ) and \ = 1 near 0, and 

D(L°°) = W 2 ' 2 (R 3 ). 
Moreover, the relationship with the previous parametrisation is 

If = L e iff ^{9) = ^ (tan(~) - l), 5 G (-7r,7r), 

and 

L°° = = A. 

Remark 2.4. Note that the case = — 7r corresponds to functions with only con- 
stant term but no ^ term in the expansion at the origin. The case 9 = | 

corresponds to functions with only ■X term but no constant term in the expansion 
at the origin. 

Proof. We focus on the first two terms in the expansion. From Proposition 12.11 
and equation (j2j) we know they are a multiple of 

(i + e *)JL_( e -?« + e <*+f)i). 

For = — 7r, the A term disappears, and we have a non-zero constant -V2. This 
means we can get any value at the origin, which is consistent with what we already 
know. In the new parametrisation, we label this operator L°°. 

Then we consider 9 ^ — n, here the first two terms in the expansion are a multiple 
of 

1 e -f i + e (Q+j)i 
\x\ l + e id 
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Denote 



A t 



(3) fi(0) = T 1 — =5 = -ei\ „ , ,. 

and in the new parametrisation this operator is denoted by ZA After some cal- 
culation, we know that 

(4) M^) = ^(tan(^)-1), 0e(-7r,7r). 

From this expression, we can see n(9) ranges across the real line. When 9 increases 
from — 7r to 7r, n{9) moves rightwards along the real line from — oo to oo, and it 
changes sign when 9 = ~. □ 

Remark 2.5. The /z-parametrization corresponds to the usual way self-adjoint ex- 
tensions of the Laplacian on a cone are defined, that is, in terms of the expansions 
of harmonic functions at the cone point, see [EMJ. 

3. Resolvent Kernel 

3.1. Resolvent Kernel. From now on we focus on M. 3 . In this section we calcu- 
late the resonances and eigenvalues of the self-adjoint extensions of L. For that 
we need to determine the resolvent kernel of these various self-adjoint extensions. 
We start with A. Since the spectrum of a self-adjoint operator is a subset of R, 
for any A ^ R, A 2 is in the resolvent set, ie (A — A 2 ) -1 exists. Let's determine the 
kernel of (A - A 2 )" 1 for A G C with Im(X) > 0. Suppose (A - \ 2 )u = f. Taking 
Fourier transform, we have (|£| 2 — A 2 )w = /. Hence, u = i|p^s , so 

(A _ xv/ = „ = ^(^) = / , r-H^) = y ^ 

Hence we know the kernel of the operator (A — A 2 ) -1 is 

e iX\x-y\ 

( 5 ) K free (X, x, y) = — r . 

4tt\x — y\ 

Now we determine the kernel of (L M — A 2 ) -1 , fx G (—00,00). Besides the term 
Kf ree (X,x,y), the kernel here has another term K extra (fi, X,x,y). We guess that 

e *A(M+[j/|) 

K extra (n, \,x,y) = b(fi, A)——-——, 

\A\y\ 

for some 6(/x, A). With this guess, we first determine what &(//, A) must be, then 
verify it is indeed the kernel what we are after. Denote 

e i\\x-y\ e i\(\x\+\y\) 

K(/j,, A, x, y) = K free (\, x, y) + K extra (n, A, x, y) = — f + &(//, A) . 

Att\x — y\ \ x \\y\ 

If K(fi, A, x, y) is indeed the kernel of (L M — A 2 ) -1 , it must lie in the domain of L M 
when y is fixed. So we fix y and consider the expansion of the function at x = 0, 

e iMy\ /ll \ 

^( 6(/i ' A) R + ^ + a6(/i ' A) ))' 
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When fi 7^ oo, b(fi, A) 7^ 0, and being in the domain of L M means 

± + i\b(n,\) 

We solve for &(//, A), 

b(ji,X) 



47r(A + i/i)' 
Therefore 



K extra {ii,\,x,y) 



i e i\{\x\+\y\) 



47r|x||y|(A + iji) ' 
and the kernel of (L M — A 2 ) -1 is 

e *A|x-j/| j iC iX{\x\+\y\) 

(6) K(n,\,x,y) — —. f + 



Att\x — y\ 4Tr\x\\y\(X + ifi) 

We are left to verify this is the correct kernel. Indeed, it maps L 2 (IR 3 ) into Z)(L M ), 
and away from the origin, we have 

(2/ - X 2 )K^ x f = (A - A 2 )i^ A / = /, 

where denotes the operator corresponds to the kernel K(/i, A, x, y). Since the 
origin has measure 0, it means in L 2 (IR 3 ) we have 

(2/ - A 2 )2^,a/ = /, 
which shows we have found the correct kernel. 

3.2. Resonance. Remember we determine the above kernel of A G C with Jm(A) > 
0. In this region, we have exponential decay, so the operator maps L 2 -functions 
to L 2 -functions. The kernel K(fi, A, x, y), as a function of A, clearly has a mero- 
moprhic continuation to the whole complex plane C. The continuation is defined 
by the same expression, so for convenience, we use the same name K([i,X,x,y) 
to denote the continuation. We see from (jBJ) that K(fi, A, x, y) has a single pole 
at A = — ifi. For fi < 0, this pole is in the physical half of the plane, and as we 
will see in the next subsection, its square is an eigenvalue of the operator ZA For 
fi > 0, the pole is in the non-physical half of the plane, but it still has physical 
significance as shown in the next section, and in this case the pole is called a 
resonance. 

Remark 3.1. Note that the Laplacian A is the only self-adjoint extension of L 
without either an eigenvalue or a resonance. 

3.3. Eigenvalue. We start with an arbitrary A, and try to find an eigenfunction 
in D(L*), then determine whether it lies in the domain of any self-adjoint extension 
of L. Suppose if is an eigenfunction in D(L*) with eigenvalue A, that means it 
lies in Ker(L* — A) = Ran(L — A) -1 , hence 

(<p, (L - A)V>) = for all ip G D(L). 

We can solve the above equation similar as before. It has a non-trivial solution 
only when A is negative, so we know none of the self-adjoint extensions has a 
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non-negative eigenvalue. For A < 0, the eigenspace corresponding to A is spanned 
by ^ 

*-'<j£pb> = sn? = m ~ tt + °<w>- 

We use equation ([2]) to obtain the expansion at the origin in the above equation. 
By Proposition [TBI we know this function is in .D(L~ v/ ~^). 

From above discussion we know that for // G (— oo,0), the extension Z/ has an 
eigenvalue — /i 2 , and the eigenspace is spanned by the function below, 

e /j.\x\ 

Since /i is negative, is in L 2 (IR 3 ). We compute that \\v^\\\ = —^7, so a nor- 
malised eigenfunction corresponding to the eigenvalue — /i 2 is 



27T \x\ 

When /i G [0, 00), that is when 6 G [f , 7r), the expression of the eigenfunction t> M 
in the above case is no longer in L 2 (IR 3 ). As mentioned earlier, in this case the 
pole a(/i) = —in is a resonance. To summarise, the spectra of the self-adjoint 
extensions are 

a(L^ = {-/i 2 }U [0,oo), /iG (-oo,0), 
er(L M ) = [0, 00), /i G [0, 00]. 



4. A WAVE EQUATION INVOLVING L M 

4.1. The wave equation. For fi G [0, 00), the resonance — ip, doesn't result in 
an eigenvalue, and it is in the non-physical half of the complex plane as it means 
exponential growth of the kernel K(fi, A, x, y). But this resonance still has physical 
significance, and as we will see, it appears in the wave kernel involving L M , see 
[LPj . The wave equation we are going to consider is 



d 2 t u + L^u = 




where /iGRU {00}, f,ge C7 6 °°(K 3 \{0}). 
We know that the solution for the system 



where a G M, is 



d\u + a 2 u = 
u\t=o = f 
d t u\ t =o = g, 



. . , sin(at) 
u(t) = cos(at)/ + — —g. 



NON-FRIED RICHS SELF-ADJOINT EXTENSIONS OF THE LAPLACIAN IN R' 



11 



So by functional calculus, the solution for the system we are interested in is 



/ r~? — \ r shifty// 1 ) 
u(t) = cos(tVL»)f + — \— - g. 

Here if / G D(L^) and g G D{\TF), we have a strong solution, ie tt(t) G D{L^) 
for each t, and it is continuous as a function of i with values in D{L^). 

We proceed to calculate the kernel of sm ^? 7 '' > , then the kernel of cos(ty/~D*) is 
given by its time derivative. 

Proposition 4.1. For am/ /iGl, i/ie kernel of hlll ^If^ for t > is 

(7) -^(5(t 2 - \x - yl 2 ) + r^n^(t - M - |y|)e^*l + '»l-*J), 

47r \x\ \y\ 

where H is the Heaviside function. 

Remark 4.2. We know that away from the origin, L M is the same as A. Also, 
due to finite propagation speed, the minimum time required to travel from x to y 
through the origin is |x| + \y\. Therefore for t < \x\ + \y\, we would expect that 

sm ^? 7 '- > has the same kernel as sin (^^) . Our kernel j7]) satifies this, and is a check 
on its correctness. 

Remark 4.3. The second term in (J7J) can be interpreted as a "diffracted" wave 
from the origin thought of as a cone point. The strength of the singularity is 1 
order weaker than the incident singularity, as is the case for a diffracted wave, see 
[CT] and [MW] . 

Remark 4.4. When /i < 0, the second term of (J7|) is exponentially growing in time 
as t — > oo, but exponentially decaying in space as \x\, \y\ —¥ oo. This is due to the 
negative eigenvalue. On the other hand, when \i > 0, this term is exponentially 
decaying in time as t — > oo, but exponentially growing in space as \x\, \y\ — > oo. 
It corresponds to a term in the "resonance expansion" for solutions to the wave 
equation on a compact set, see |LP] and |TZ] . 

Remark 4.5. See |LH] for a different approach to obtain the kernel by solving an 
auxiliary problem. 

Proof. We have, 

Mt^)_ rw[ dP ^ lim r-n^^ 



where (p : R — > [0, 1] is a smooth function such that (p = 1 on -&l(0) and (p — 
outside -62(0), and the limit converges under the strong operator topology. Then 
depending on whether Z/ has an eigenvalue, we have two possibilities. First for 
H G [0, 00), we have 



*M = lim rMtv^ ^ lim r-E&^w. 
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While for oo, 0) the eigenvalue contributes an extra term, 



sin(tVi>) ,. f°° sm(ty/a) ,<r Sjrt .. r°°wa.(ty/F\ ,<r. Jrt sin(i/xt) „ 
-lun / — )gLP ct = hm / — v Jl )c?P (T H y~-Lp_ 



2, 



where P-^ is the orthogonal projection onto the eigenspace of — /i 2 . 

In either case we need to calculate the term lim^oo J °° sm ^^ (p(^)dP a . For each 
R > 0, we apply integration by parts twice to evaluate the integral, and also by 
using Stone's formula, we obtain, 



Sin( ^(-)^ 



1 f°°sm(ty/a) ,a. v 



r ™P*>(~) Hm ((I/ 1 - a - ^ - (I/* — a + ie)^)da. 
Jo v a e +° 



2vri 

Remember the kernel of (D* — A 2 )" 1 , // G R U {oo} is dented by K(fi, X, x, y), 



hence for any function / G D{\/L^) we have 



I sin (t o- ^o. I ^ K (j 1 ^^ x ^_x{fi i -y/^ 7 x,y))f{y)dycb 

JO V a n J-oo 

roo poo ^2 

— / / sin(tA)iT(/i, A, x, y)tp( — )f(y)dXdy (substitute A = y/a). 

7n y_oo y^oo R 



From before, we know that 

K(n,\,x,y) = K free (X,x,y)+K extra (ii,X,x,y) = ^ - ^ + ^\ x \\yU^ + ifi y 

The free resolvent kernel Kf ree (X,x,y) substituted to the last line of (|SJ) gives the 
free wave kernel -^S(t 2 — \x — y\ 2 ), so from now on we concentrate on the term 
contributed by K extra (fi, A, x, y), 

1 f°° A 2 

(9) — lim / sm(tX)K extra (n,X,x,y)(p(—)dX. 

TCI R^-oo J _ QO K 

We continue the computation in cases depending on the sign of \i. 

Case 1: when \l — 0. 

In this case equation (jUJ) becomes 

i r 00 1 \ 2 



47T 2 |x||?/| R- 



/oo i \2 
-sin(tA)e iA(|a ' l+|y| V(— )dX 
oo A -R 



(10) 



1 f°° 1 A 

, a, I, | J im / T sin(tA)cos(A(|x| + \y\))<p(— )dX 
Arc \x\\y\ R-^ 00 J -oo A ' R 

i / r°° i a 2 

— n . n i hm f / — sin (X(t + \x\ + \y\))ip(— )dX 

87T 2 |x||?/| R-^oo V J_ 00 A V V 11 W'^rJ 

/°° 1 A 2 
-sin (A(t-|x|-|y|))^(— )dX 
-OO A ^ 
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We then make a substitution, and split into three cases: 
(i) When t — \x\ — \y\ > 0, ( |T0l) becomes 

(ii) v 

note that we have two Dirichlet integrals, each of which equals 7r, so the 
above expression equals 

1 



Atc\x\ \y\. 

(ii) When t — \x\ — \y\ < 0, after the substitution we have the same expression 
as (TTTT) except the sign of the second integral is negative. Hence the two 
integrals cancel each other, therefore in this case (llOp equals 0. 
(hi) When t —\x\ — \y\ = 0, the second integral in the last line of equation ffTUl) 
is 0, hence we only get the first integral in expression (TTTT) . therefore here 
(TO) equals 

Combine all three cases, the kernel of sm ^0~^ contributed by K extra for t > is 



' -Mt-\x\-\y\). 



An\x\ \y\ 
Case 2: when [i ^ 0. 

Since sin(tA) = e ' ~^ , we have to deal with the following two terms 

r J.\{\x\+\y\+t) \2 

(12) „, ,, , lim / <p(—)d\, 

8ir 2 \x\\y\ R->°° J Im (\)=o X + i/j R 

and 

• r i\{\x\ + \y\-t) \2 

(13) „, ,, , lim / f(^)dX. 

%ir 2 \x\\y\ R-+°° J Irn ( X)=Q X + i/j R 



We first deal with the integral ((T2 



r e iX(\x\+\y\+t) ^2 

lim / — ; ip(—)d\ 

R ~+°° Jlm(X)=0 A + Ifl R 



(14) = - , \ , lim / A( e iA(l»l+l»l+t))^4ldA 



i{\x\ + |y| +t) ^i; ra (A)=o dA v " 7 A + z// 

/ A 2 \ 



i(\x\ + \y\ + t) R-^c 



lim [ e ^\ + \y\ + t)±^^Rl) dX . 
R ^°°Ji m (X)=o dX X + ifi' 

The last equality is established by integration by parts. Apply the quotient rule 
then we get two integrals. The first one is 



2 f \ v > ^yH\*\+\y\+t) 



lim — / H ' dX. 

R-^oc R J Im{x)=Q X + ifi 
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This limit is 0, hence ( I14p equals the second integral obtained from the application 
of quotient rule, which is 

lim / ; - R dX 



1 


i{\x\ + \y\ - 


ft) 


1 




i(\x\ + \y\ - 


ft) 


1 




i(\x\ + \y\ - 


ft) 



/m(A)=0 (A + lUf 

e^H+lwl+Olim^^) 

^ d ^ 

:0 (A + iii) 1 

(by Dominated Convergence Theorem) 

e i\(\x\ + \y\+t) 



-dX. 



Therefore term (|12|) becomes 

_1 r e i\(\x\+\y\+t) 

(15) 87r 2 |x|M(|x| + M+t)/ /m(A)=0 (A + z/i) 2 dX ' 

Now we can shift the contour Jm(A) = upwards to Im(X) = M for any M > 0, 
and the integral should stay the same except when the contour moves across a 
pole. When M — > oo, the integral approaches zero. 

Similarly, term ffTS]) becomes 

1 /• e iA(|s|+|tf|-t) 

(l(>) 87r 2 |x||y|(|x| + \y\ -t) J Im{X )=o (A + i/i) 2 dX ' 

In the region t < |x| + we can shift the contour Im(X) = upwards, while 
in the region t > \x\ + we can shift it downwards. As before, the integral 
stays the same except when the contour moves across a pole, and the integral 
approaches zero when the contour is shifted further and further away. Note that 
for t = | x | + \y\, we have a choice between shifting it upwards or downwards so 
we can always avoid the pole, hence we know the integral is 0. We continue the 
computation in two subcases: 

Subcase 2(a): when fi > 0. 

In this case the pole, which is the resonance a(/i) = —in, lies on the negative 
imaginary axis. We shift the contour upwards for the integral (I15p . and the integral 
goes to 0. While for the integral ( ITBjl . it depends on the sign of |x| + \y\ — t: 

(i) In the region t < \x\ + \y\, the contour is also shifted upwards, so it also 
goes to 0. Hence equals 0. 

(ii) In the region t > \x\ + \y\, the contour is shifted downwards hence across 
the pole —ill. The residue of the integrand at the pole A = —ift is 

i(\x\ + |y|-t)e^ l+l2/| -*), 

since the winding number is —1, by Residue Theorem, we know f lT6|) . ie 
([9]) equals 

e MM+li/M) 
4n\x\\y\ 
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Combine the two cases using a single expression, the kernel of sm ^f^ contributed 
by K extra for \i > and t > is 



H(t- 


x\ 




\y\)e 


A\*\+\v\-t) 


47T 


x\ 


\y\ 





Subcase 2(b): when ji < 0. 

In this case we have a negative eigenvalue and, 

sin(tv^) Um f°° Mtf)^ )dPa + BinM) p _ 



Let's first deal with the E£i!H*2 p 2 term. To do that we calculate the kernel 
of P_ M 2. As discussed before, the eigenspace of — /i 2 is one dimensional, and a 
normalised eigenfunction is 



v^x) 



2tt \x\ 



so the kernel of P_ M 2 is 



_^ e MM+M) 
2vr|x| \y\ 



v^x)v^y) ■■ 

Therefore the kernel of ""l 1 ^ P_.,a is 

sin(z//i) -^e^M+lvl) _ sinh(^t) -^eMkl+lwl) _ _ e A»(M+[vl+t) ^ e M(M+[vM) 
z/i 27r|rr||y| // 27r|x||?/| 47r | a: | 47r|x||?/| 

Now we calculate lim^oo J °° cos(ty/a)(p(^)dP (T . The pole, ie the resonance a(/i) = 
—in, now lies on the positive imaginary axis, hence integral ( fl5|) may contribute 
some value. The residue of the integrand at the pole A = — i\x is 

%{\x\ + \y\ +t) e "(M+M+*), 
since the winding number is 1, by the Residue Theorem, the integral ffl5|) equals 

( 18 ) 

47r|s| \y\ 

Again as before, to deal with ( fl~6|) we split into two cases. 

(i) In the region t < \x\ + \y\, the contour is shifted upwards, hence across the 
pole — ifM. The residue of the integrand at the pole A = —i\x is 

i(\x\ + \y\ -t) e "(H+W-*). 
The winding number is 1, so by Residue Theorem, we know ffl6l) equals 

- e ti(M+\y\-t) 



(19) 



4tt\x\ \y\ 

Adding (1T71), (fig]) and (1T9]) gives us 0. 



16 



PAUL LIN 



(ii) In the region t > \x\ + \y\, the contour is shifted downwards, so it doesn't 
move across any pole, so ffTB"]) is here. We add (IT7|) and ffl~8]) to get 

el i(\x\ + \y\-t) 

4vr|x| 1 2/ 1 

Combine above two cases into a single expression, the kernel of sm ^r^ con- 
tributed by K extra for /i < and t > is 



H(t- 


x\ 




\y 




A\x\+\v\-t) 


47T 


x\ 


\y\ 





We have completed the proof. 

□ 

Remark 4.6. In the three cases of different signs of /i, we have different expressions 
for the kernel for t = \x\ + \y\. This is fine because we can ignore the value of the 
kernel on a set of measure 0, but we do need to know that it is finite there so that 
we can be sure there is no distribution supported on this set. 
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